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MODULATION AND SAMPLING OF HYDROMAGNET IC RADIATION* 


By Charles P. Sonett 
Ames Research Center 

SUMMARY 

It is shown that the detection of electromagnetic signals in 
the hydromagnetic domain by means of magnetometers is a low fre- 
quency extension of the electromagnetic spectrum; therefore con- 
ventional modulation theory is employable , subject to the special 
constraint that the carrier frequency is determined by spacecraft 
motions. In the presence of a hydromagnetic continuum, spacecraft 
spin modulation folds the spectrum so that fictitious offsets , 
ambient field errors, and the generation of spurious wave signals 
can occur, the latter from single sidebanding. Synchronous demod- 
ulation by means of time quadrature detectors can be employed to 
remove such effects . 

Commutated telemetry which time shares experiments is shown 
to yield a multirate sampled spectrum. Thus telemetry can cause 
uncorrectable spectral folding of many orders because it generates 
a complex switching spectrum containing a subsidiary multiplet 
structure about each of the principal switching tones. 
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1 . Intr oduct ion 


The magneto spheric and interplanetary plasmas are regarded 
as fluids to which is coupled an ambient magnetic field with an 
electromagnetic radiation field completing the ensemble . The 
spectral content is determined locally both by plasma kinetics 
and the form of the source which drives the time varying field, 
except when local plasma instabilities are important, or the 
decay of turbulent motions is dominant. The rapid flow of the 
interplanetary plasma against a spacecraft means that all fre- 
quencies are Doppler shifted and time versus space variations 
cannot be simply distinguished. Thus the time reversibility of 
propagating waves becomes mixed with the diffusion-like and 
irreversible decay of large spatial disturbances into smaller 
ones, the familiar cascading from small to large wave numbers 
in turbulence . 

Of special interest is the low frequency region of the 
spectrum associated with hydromagnetic radiation, as both wave 
and convective frequencies of the solar gas encompass character- 
istic frequencies of motion of spacecraft. It ensues that cer- 
tain modulation effects take place which can modify the 
interpretation of signals from the field. The modulation has a 
geometric basis which in principle generalizes to any detector 
in motion. Spacecraft rotation has an especially simple modu- 
lation effect upon magnetic signals which can be explained by 
an extension of existing modulation theory. The special quality 
of the modulation arising from spacecraft spin is due to the 



superposition of similar frequencies from the radiation field 
and the spacecraft. This situation favors the development of 
asymmetric and truncated side hand spectra which in turn can 
"alias" the signal.* Variations of the discussion of this 
paper can he extended, for example, to particle detectors 
following the arguments laid down here. 

The plan of the paper is to consider first what happens 
when a magnetic field in space is measured hy a magnetometer 
mounted upon a spacecraft which is executing motions about its 
center of gravity. Specifically, the magnetometer will he 
assumed to decompose the field into three orthogonal components 
as is the case for the triaxial flux gate. The spacecraft will 
he assumed to spin about the figure axis, as this analysis pro- 
vides a realistic view of the resultant modulation process. 

Both the changes in amplitude and changes in direction of the 
field are examined separately, and the general case of both 
types of time dependence occurring simultaneously will he shown 
to he a simple combination of these two. For instruments which 
measure the field directly, in contrast, say, to search coils 
which operate upon the field hy time differentiation, the two 

*The term alias is jargon used to denote first order fold- 
ing of spectra due to violation of the Nyquist limit in the 
sampling of continuous data. This special case, discussed later 
in the text (see fig. 11 ), is used interchangeably in this paper 
with spectral shifting which is a generalization of first-order 


spectral folding . 



distinct types of field variation generate, respectively, 
amplitude and phase modulation of the output signal from the 
spacecraft. For a nonrotating spacecraft platform, the former 
degenerates to the trivial case of AM modulation with zero 
carrier frequency whereas directional changes in the field gen- 
erate more complex signals. For brevity these two effects are 
termed scalar and vector modulation, respectively, recognizing 
that the latter term is not exact. This discussion is followed 
by a discussion of the search coil on a spinning spacecraft. 

The signal properties are shown to be different according to the 
operation of time differentiation inherent in this type of 
detector. Following this, attention is directed to the corre- 
spondence of data multiplexing (customarily employed in single 
channel spacecraft telemetry) to multirate switching and the 
spectral overlap or alias problem. Lastly, the application of 
synchronous demodulation to spin modulated signals is considered. 
Synchronous detection employing quadrature signals provides a 
means for the electrical despinning essential for the unfolding 
of side band structure. 

In the discussion that follows it is always assumed that 
the system response is linear, so that harmonic generation from 
this source is absent. Thus the arguments entail only modulation 
of a linear system, and spectral folding and shifting from multi- 
periodic sampling as occurs in time shared telemetry. Also argu- 
ments regarding spacecraft motions which involve precession and 
nutation are not considered, having been discussed elsewhere [1,2]. 
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2.0 Convection and Radiation 


It is thought that in the solar wind, spatial structures 
of nearly any size should he allowed and the gyro radius is not 
an essential limit, though the more time stationary and organized 
forms such as shocks and contact discontinuities should display 
a limit more or less corresponding to the ion gyro radius. Since 
the wind is in motion with a proper velocity of some hundreds of 
kilometers per second, it is clear that such a limit means that 
to observe the details of diffusion of a contact discontinuity 
or other similar organized change in the gas and the field requires 
that a magnetometer have a frequency response of at least 1 cps 
or more • There are other phenomena in the interplanetary gas which 
are presently little understood such as the alternating polarity 
of tubes of force emanating from the sun. How these merge where 
opposed polarities meet has at least superficial similarity to 
questions relating to contact discontinuities. 

Although propagating waves are, in practice today, indis- 
tinguishable from spatial changes in the solar wind, in the mag- 
netosphere bulk velocities should be much reduced and perhaps 
absent from time to time, though there is no empirical basis for 
assuming this. In either case when the field is sufficiently 
steady and large to form a reference frame for waves, the geom- 
etry rules governing plasma waves should be applicable to con- 
vected disturbances. 

Only a simplified view of wave phenomena is required for the 
ensuing arguments as they will rest partially upon geometry. 
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Consideration of the allowed wave modes is restricted to the 
hydromagnetic domain governed primarily by ion motions. Exten- 
sion upwards implies inclusion of frequencies not properly the 
province of hydromagnetics. From geometrical considerations 
alone the three governing modes become the Alfven shear wave 
and the two mixed waves representing fast and slow modes. These 
have phase velocities which depend both upon the Alfven and sonic 
speeds. The corresponding surfaces of constant phase are shown 
in fig. 1. It is seen that propagation across the field is 
allowed only for the fast mode while along the field direction 
two waves can propagate. The diagram is the zero frequency 
limit and displays no dispersion. (in a broader sense , however , 
there is a functional relation between speed of propagation and 
direction with respect to the preferred axis of the field.) 

The linearly polarized Alfven mode displays rotation of the 
plane of polarization as the gyro frequency of the plasma is 
approached from the lower side. Above u>g only one -handedness 
persists. Thus the Alfven wave is circularly polarized in 
regions near to and above the gyro resonance . In interplane- 
tary space tOg is typically of the order of 1 cps and therefore 
is important for magnetometer experiments . For the Alfven mode 
the disturbance vector is normal to the field and gyrates about 
the latter, as shown in fig. 2. Propagation across the field 
consists of alternate collective compression and rarefaction 
in the field and the gas pressure, while for the Alfven mode 



the motion is incompressible. These two limiting conditions 
form the geometrical background for the discussion of modulation 
effects to follow. 

Purely compressive changes in the field associated with 
progressive or standing waves are subject to amplitude modula- 
tion by the spinning reference system. Alfven waves are phase 
modulated by the spin. It is more convenient to invert these 
problems so as to regard the spin as the carrier and the hydro- 
magnetic radiation as the source of the modulation. There is 
no loss of generality in doing this and the analysis, in par- 
ticular, of the phase modulation is simplified considerably. 

By way of summary, it is assumed that variations in the 
field, both spatial and time dependent, are representable by 
simple geometric constructions so that the main points discussed 
are not obscured by detailed development of the allowed modes 
for propagating waves in a plasma. The problem becomes one of 
synthesizing the signal characteristics of magnetic field vari- 
ations from data taken in a hydromagnet ic radiation field hav- 
ing spectral components which beat with the spin frequency of 
the spacecraft. Since spatial variations are not separable 
without recourse to other data, especially plasma probe infor- 
mation, for the purposes of this discussion they will be treated 
as strictly time variant. Further we shall assume that the 
hydromagnetic field consists of a continuum of frequencies, and 
only for the purposes of illustration will deviation be made 
from this choice . 
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The compressive case corresponds to the rms sum of the 
Alfven and the sonic speeds taken across the preferred direc- 
tion of the field. This is the most straightforward case to 
treat as the field suffers only changes in amplitude and thus 
the magnetic signal consists of heats between the spin and the 
hydromagnetic spectrum. The preferred direction of the mag- 
netic field is generally not along the spin axis of the space- 
craft; thus the field is resolved into spacecraft coordinates. 

3-0 The Scalar Hydromagn etic Radiation Field and Spectral Folding 

The simplest radiation field to treat is that corresponding 
to the magnetic field undergoing compressive changes associated 
with fast wave propagation across the field lines. This simple 
time variance is now considered with the signal detected upon a 
spinning spacecraft, and it will be shown that a spectral line 
splits in a manner strictly analogous to amplitude modulation 
where the spin frequency of the spacecraft assumes the role of 
the carrier. 

Since the field static direction is assumed known, as is 
that of the spacecraft , the components of the field are resolved 
into spacecraft coordinates directly. Let (xi,x 2 ,x£) be a right- 
handed, mutually orthogonal set of unit vectors associated with 
the ambient field passing through the spacecraft position with 
Xi along the line of force taken positive in the direction of 
the field, x 2 the unit vector lying in the magnetic meridian 
plane and having positive sense outward (for the interplanetary 
field where a preferred direction is more difficult to define 



other appropriate coordinates such as solar ecliptic are 
available), and X3 the -unit binomial to the field line. The 
notation (xi,x 2 ,Xq) is given to a second set of unit vectors 
centered in the spacecraft but nonrotating. The spin axis of 
the vehicle is given by w s = Jw s |x3 which defines the direc- 
tion of the unit vector , X3 . It remains to define the line of 
nodes between the two sets of vectors to complete the descrip- 
tion.* The orientation of the spacecraft with respect to the 
magnetic field is then given by the two Eulerian angles, G and 
cp. The third angle, 3|r, the rotation about the axis, X3, is 
arbitrary and taken to be zero. The nominal or ambient mag- 
netic field is given by H 0 = |H 0 |xi. To this must be added 
the perturbation field which is represented by the Fourier 
transform pair 

H(t) = (l/y/Et) p h(w)exp(iwt)dw (l) 

00 


H(w) = (l/V2rt) P° H(t)exp(-icot)dt (2) 

J -00 


where H(t) is the disturbance waveform and h(u>) is the spectral 
density of the field disturbance in units of field intensity 
per radian per second. For the case at hand where the field 
variations are taken to be strictly compressive, the variation 
field is in the direction Xi. Thus the total field is given by 


Hi = H 0 + H(t) 
H 2 = 0 > 

H 3 = 0 


( 3 ) 


*The Euler angle definitions used follow those of Goldstein [3 3 • 



In the nonrotating , spacecraft centered coordinate system 
(xi ; x 2 42 ); the components of the field are resolved into 

a Hi = Hi cos 9 + H 2 sin cp (4a) 

a H 2 = -Hi sin cp cos 0 + H 2 cos cp cos 0 + H 3 sin 0 (4b) 

a H 3 = Hi sin 9 sin 0 - H 2 cos 9 sin 0 + H 3 cos 0 (4c) 

The angles 0 and 9 are taken to be fixed for times comparable 
to those of interest . Thus in turn ^H 2 and ^E 3 may be regarded 
as fixed. 

3*1 The Generation of AM Side Bands 

Consider now a triaxial set of sensors, such as flux gates 
mounted upon a spacecraft having spin angular velocity w s . The 
sensors lie, respectively, along the spin axis, S3, and mutually 
orthogonal in the equatorial plane. Let 61 and c 2 be the unit 
vectors associated with the two equatorial sensors. The time 
origin is chosen so that c = ci cos(co s t) + c 2 sin(w s t) . Then 
the components of the field as observed from the rotating plat- 
form of the spacecraft are given by (H • c) or 

^Hi cos(w s t) ( 5 a) 

A H 2 sin(w s t) (5b) 

A 

where the phase is arbitrarily adjusted so H = a Hi at t = 0. 

For algebraic simplicity, the spin tone is collapsed to a 
delta function and convolution of the spin tone with the hydro- 
magnetic field yields a total spectral density distribution 
given by 


H q 5(u) s ) + (l/ 2 )h(w) [6(u) b + w) + b(w s - w) ] 


(6) 



which means that the spin tone heats with the hydromagnetic 
spectrum corresponding to direct amplitude modulation- The 
special quality is associated with the embedding of the spin 
tone or carrier frequency in the frequency domain of the hydro- 
magnetic radiation field. The second two terms in eq. (6) 
together with the carrier can he represented in the conventional 
manner as spinning vectors as shown in fig- 3* In the inertial 
frame the two vectors rotate with angular velocities co s + u) 
and u) s - w. In the frame which transforms away H 0 , they rotate 
with velocities u) and -u>. This is an important effect for 
the situation discussed in the next section. 

Contributions to h(oo) neighboring u> s are translated to 
regions near the origin while those near the origin translate 
to the neighborhood of w s . Alternatively on a frequency plot 
showing negative frequencies, side band symmetry about co s 
would extend the spectrum to negative frequencies. The mixing 
of the hydromagnetic spectrum with the spin tone splits the power 
shown in the side band structure according to eq. (6) where half 
the power from a band is shifted into each of the side band 

components. The overlapping of translated side bands, however, 
does not allow a direct computation of the power in a particular 
region of the spectrum unless the functional form of h(w) is 
given . 

When the frequency band lies in the interval 0 < Aa) < u) s 
the upper side band is found at w > w s . For frequencies in the 


interval w s < Aw < 2to s , reflection about zero frequency occurs 
so that the lower side band is still contained in the interval 
0 < Ah < w s with the spectral form inverted. However, for 
those regions of the frequency domain of h(w) which lie above 
2w s , the modulation process generates both side bands asymmetri- 
cally above the carrier. This is a direct result of the low 
carrier frequency and is never found in normal communication 
practice [ 4 ] where the modulation is bounded to a small frac- 
tion of the carrier frequency (see fig. 4 ). 

3 *2 Single Side Band Ge neration of Pseudo-Alfven Signals 

Consider a real system where the upper frequency limit is 
defined by a filter that is flat out to co c where w c > w s 
and rolls off beyond w c at some preassigned rate. Then al 1 
frequencies below w* < w c - u> s will be recovered without 
spurious information. For those contributions to h(w) where 
w* > to c - u) s , the upper and lower side bands will display 
different amplitudes according to the expression 

H 0 5 (u) s ) + (l/2)h(to) [a6(u) s - w) + b6(w s + w) ] (7) 

where a is unity for modulation frequencies, u> < u s + co Q , 
h is unity for w < w c - u) s , and h(w) is the spectral density 
when w < u> c . When h(w) lies in the range of the slope of the 
filter, then it must also he replaced hy an equivalent h*(w) 
since the modulation tone is attenuated. From eq. (7) it is 
seen that the side hand components are attenuated hy different 
amounts since they lie at different points on the filter ramp. 
Reference to fig. 5 a shows an idealization of the above. 



An important effect of the truncation of the side hand 
structure is to introduce a spurious phase modulation whose 
most general interpretation suggests the presence of Alfven 
waves when none are present in the field. This is seen 
hy reference to fig. 5h where a vectorial representation of 
the pertinent signals shows the manner of composition of the 
total signal. When the side hand structure is symmetric, it 
is seen that the two side hands are symmetrically disposed 
about the carrier vector, with the result that phase contri- 
butions are cancelled. 

The case of the line spectrum sho™. in fig . 5 has the 
equivalent time expression (see fig. 6) 

H Q cos u) s t + (cos w s t ■ cos w-t)+ (sin w s t • sin 

( 8 ) 

and the quadrature side hands are in ratio (a + b)/(a - h) . 

When the upper side hand is suppressed completely, the in- 
phase and quadrature components are equal and the power in 
the sinusoidal and cosinusoidal parts of the carrier are 
equal. This means that only in the event the value of h(w) 
drops to zero for all a) > w c - oo s can it he guaranteed that 
a spurious Alfven wave field will not he introduced. The 
results of this section can he summarized as follows: 

In the presence of a scalar hydromagnetic radiation field , 
that is, containing only waves characterized hy changes in 


field amplitude, t he signal from a magnetometer sensor which 



measures the instantaneo us field normal to the spin axis of a 
rotating spacecraft h aving a spin frequency, gj s , will contain 
loth in-phase and q uadrature modulation side bands derived from 
that portion of the ra diation spectrum lying at frequencies in 
the range io c - ws < 10 < Wq + wg where w c is the filter break 
frequency of the system . 

3 • 3 Spacecraft Pseudo-Per m and Amb lent Field Offset 

Consider now that portion of the hydromagnetic spectrum, 
h(w), that overlaps the spin frequency, w s . Because of modula- 
tion, this portion of the spectrum is folded to the neighborhood 
of to = 0; consequently a pseudo-bias appears in the spin signal 
as perceived at the receiver. Depending upon the time of averag- 
ing, which determines the line width at zero frequency, an error 
could be introduced into the measurement of the offset which, 
in turn, measures the component of the spacecraft bias lying in 
the equatorial plane of the spacecraft. Although in practice 
this offset of the magnetometer signal may be undetectable, it 
is important to point out that in principal it is always present. 
Its value will depend, as stated above, upon the filter band- 
width used to assess the offset and the power in the hydromag- 
netic spectrum about the spin frequency. 

The arguments just given can be applied in a broader context. 
For example, if the hydromagnetic field contains frequencies 
which are dominant at twice the spin frequency, the folding of 
these with the spin frequency will cause them to appear as addi- 
tions to the basic sinusoidal signal. Customarily this signal 


is interpreted as arising solely from the static component of 
the magnetic field which is projected into the equatorial plane 
of the spacecraft. Then the measurement of the static field will 
he increased hy some amount determined hy the power in the hydro- 
magnetic domain nearly twice the spin frequency and the length 
of the sampling interval used to measure the spin tone, or, 
equivalently, the filter bandwidth. 

In arguing the importance of these folded contributions to 
the measurement of the static field and the spacecraft perm, it 
is clearly important to understand the form of the spectral com- 
ponents involved in the process. To a point, common sense must 
be applied in studying data of the type discussed. However, 
cases can be expected to arise where, in spite of the natural 
appearance of the data, it does happen that important folding 
from the modulation takes place; such cases cannot be distin- 
guished from the unfolded data by any a posteriori means . 

3.4 The A periodic Field Reversal 

Occasions arise especially in the interplanetary field 
when a reversal in polarity takes place over a time interval of 
the order of five minutes by present best estimate [5^6]. 

Actually the time interval may be considerably shorter; it is 
also important to recognize the role of convection in this 
description since the field Is being swept past the magnetom- 
eter. Thus the change of field may be primarily spatial, but 
in the coordinate system of the magnetometer, can be considered 
to be strictly a time variation. Interest stems as before in 
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the role played hy folding of the spectral components of the 
field change with the spin tone from a rotating sensor. A 
quantitative treatment is beyond the scope of this paper; 
however, the simple field discontinuity can be represented 
by a step function; it is recognized that more complicated 
field changes can take place involving skewing, with the field 
never reaching zero magnitude. Such changes are also outside 
the province of this discussion as they involve rotation of 
the magnetic vector corresponding to phase modulation discussed 
in the next section. 

The Fourier transform of the step function is of the form 
l/oj and is unbounded near the origin. Thus when folded onto 
the spin frequency, the net spectrum takes on the character 
displayed in fig. 7 which shows the low frequencies shifted to 
the spin tone. The regions near zero frequency are not espe- 
cially interesting from a physical standpoint as they correspond 
to the tails of the function where the value is steady. Of 
principal interest are the higher frequencies which provide 
information on the thinness of the discontinuity and on fluc- 
tuations. In spite of the divergence of the transform around 
the spin frequency as shown in fig. 7; the power must corre- 
spond to the steady-state field far from the discontinuity. 

There is considerable spectral energy near zero frequency which, 
correctly interpreted, should correspond to regions of the step 
spectrum near the spin frequency. Lastly, there is a strong 
increment of the spectrum which is shifted to higher frequencies. 



The general effect of these folding phenomena is to decrease 
the apparent thickness of the region of transition, to possibly 
introduce a pulse-like change in the spin tone and, lastly, to 
introduce an aperiodic transient change of the zero or average 
level of the magnetometer signal. 

4.0 Vector Modulation 

In this section the signal resulting from the detection of 
the propagating Alfven mode is considered. The geometry of the 
wave in the previous coordinate system (xi,x 2 ,Xs) will be the 
basic reference for the following discussion. We assume that 
the mode is linearly polarized but that the plane of polariza- 
tion turns slowly subject to an ion Faraday rotation. Wear 
the gyro resonance for protons, one of the modes is extinguished 
and the wave is regarded as circularly polarized with a gyration 
period equal to the Larmor period. Since the ambient field 
has a fixed orientation, the resolution of the wave into com- 
ponents of the prescribed set follows directly the rules for 
phase modulation with suppressed carrier. Thus the field com- 
ponents are given by 


Hi = H Q cos | sin(& sin vt) 
H 2 = H 0 sin £ sin(& sin vt) 
H3 = H Q cos(s sin vt) 


(9) 


where | is the dihedral angle "between the instantaneous plane 
containing the wave and the plane X1X3. The other quantities 
are defined as follows: v is the frequency of the Alfven wave 



assumed monochromatic, and 5 is the deviation in phase 
corresponding to the angular deviation of the quiescent field 
around its nominal orientation. From eqs. (9) it is clear 
that the geometrical resolution of the time-dependent field 
directly generates phase modulated signals in each of the 
sensors. In general, the spacecraft will not have a preferred 
direction and therefore the components of eqs, (9) must he 
resolved through an Eulerian transformation. In the new set 
(xi,x 2 ,%3) defined as in the previous section centered in the 
spacecraft, the field is given hy 

H = [T]H ( 10 ) 

where [T] is the Eulerian rotation matrix. In the most gen- 
eral instance the spacecraft is allowed to rotate about its 
axis of symmetry with = co s t. More general motions described 
as tumbling where all the Eulerian angles become time dependent 
are not considered. In the instance of spin the process of 
signal modulation enters as in the earlier discussion of ampli- 
tude effects though in a more complex way. Even the resolution 
into components itself is important in mixing the modulation 
components. The components of the field under phase modulation 
of deviation, 6, at frequency, v, are expressed for a general 
orientation by 

a Hi = H 0 [an cos | sin (5 sin vt) 

+ a 22 sin £ sin (5 sin vt) + a 13 cos(& sin vt)] (lla) 
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A H 2 = H Q [a 21 cos | sin(S sin vt) 

+ a 22 sin £ sin(S sin vt) + a 23 cos(S sin vt)] (lit) 

A H 3 = H Q [a 3 i cos £ sin(& sin vt) 

-i- a 32 sin £ sin(5 sin vt) + a 33 cos(5 sin vt)] (lie) 

where a Hi, A H 2 , and A H 3 are the components of the field in the 
~ system, and the aqj are the elements of the matrix, given in 
eq. (10) . 

It is to he noted that linear mixing of the various side hand 
tone systems depends on hoth the orientation of the plane of polar- 
ization of the wave and the orientation of the spacecraft as 
reflected in the elements of the transformation matrix, a*^. 

The terms in eqs . (ll) are recognized as the usual expressions 
encountered in phase modulation. The carrier is suppressed and 
a Hi, a H 2 , and A H 3 contribute combinations of odd and even har- 
monics. The amplitudes of the side hand tones are given by 
Bessel functions of the first kind, utilizing the following 
expansions: 


00 

sin(S sin vt) = 2^ J^ +i (&) sin(2n + l) vt (12a) 

n=-oo 

00 

cos(5 sin vt) = J Q (5) + 2 ^)T J 2n^ cos (12b) 

n=- oo 

When these are inserted in eqs. (ll) it is recognized that for the 
most general orientation of the spacecraft, each of the component 
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magnetometers measures a spectral complex consisting of linear 
combinations of odd and even terms with corresponding Bessel 
function amplitudes. The spectrum of each of the expressions 
given above is shown in fig. 8 for the case chosen of a maximum 
phase deviation of 1.5 radians corresponding to a large Alfven 
wave which causes an instantaneous variation in field direction 
of Also it is assumed in this figure that a wave train 

exists with a single frequency given by v. The number of sig- 
nificant side bands depends only upon the phase deviation and 
therefore the total bandwidth encompassed varies directly with 
the phase excursion frequency; that is, the separation between 
side bands varies with v. Since the modulation is about zero 
frequency, i.e., the static field, the field aligned spectral 
components are even functions; fig. 8 represents terms which 
when folded add to double value in the positive frequency 
domain. When the power spectra are considered, only the posi- 
tive frequency domain is present. The components from the 
negative domain are folded about zero frequency. Mathematically 
this follows from the odd and even order Bessel functions being 
of odd and even symmetry, respectively. Thus in eqs . (12a) 
and (12b) the expressions are seen to be even for each case. 

*In all expansions discussed in this section the frequency 
domain is taken over -oo < u < oo for reasons of symmetry. The 
practical case of measurement discussed in section 7 reverts to 
folding about w = 0 since laboratory frequency analysis precludes 
a) < 0. 
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In the event that sufficient power is available in the 
side band structure for detailed analysis it becomes evident 
that field orientation information can be derived from the 
relative amplitudes of the even and odd side bands. These 
determine the orientation of the figure axis of the spacecraft 
with the direction of the ambient field. Regarding the sup- 
pressed carrier condition, it is seen by examination of 
eq. (12a) that since the signal spectrum is expanded about zero 
frequency, the carrier is nonexistent for the odd harmonic 
series. However, carrier suppression for the even series can 
also occur when the field is deviated by 2.4 radians which is 
the first root for J Q . Zeros for the higher order functions 
would be physically unreal and are therefore not discussed. 
Even deviations of 2.4 radians seem unlikely for a physically 
real Alfven wave as the field would have to be folded over on 
itself.* We now write the complete expressions for the field 
components for the three axes for arbitrary orientation of 

the spacecraft. Erom eqs . (ll) and (12) they are 

f r oo i 




' — — 


J 2n +i^ sln ( 2ri + !) vt 


n=-oo 


+ a i3 


Jo(8) + 2 


J PT1 (S)cos 2nvt 


n — oo 


(13a) 


*rr 


The restriction in angular deviation would not necessarily 
hold for the case where aperiodic changes in direction of the 
field occurs as where tubes of opposed polarity swept across the 


spacecraft . 



,H 2 = H 0 


(a 2i cos g + a 22 sin |)x 2 ) J 2n+i (5) sin(2n +l) vt 


n=-oo 


+ a 23 


,(5) + 2^ J 2n (6) cos 2nvt 


n=-oo 


(13b) 


,H 3 = H r 


00 

(a 3 i cos g +a 32 sin g ) X 2^' J 2n+1 (5) sin(2n +l) vt 


n=-oo 


+ a 33 


-00 — 

J Q (6) + 2^ J 2n (S)cos 2nvt 


n=- oo 


(13c) 


These are the most general expressions for pure Alfven wave 
modulation. Below about 0.3 radian deviation the side band 
structure collapses to two lines given by the terms Jj.. Thus 
for small amplitude Alfven waves the modulation spectrum is 
barely distinguishable from that produced by simple scalar modu- 
lation of the field. The separation of the side bands is still 
given by the modulating frequency. The even order series of 
terms vanishes and the scheme becomes especially simple. The 
only feature distinguishable from the scalar modulation is _that 
the upper and lower side bands are in phase opposition, and 
modulation repre sented by the even series vanishes. 

It is anticipated from the foregoing discussion that the 
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internal bandwidth requirement for nontruncated detection 
of vector modulation can be greater than for the scalar 
case when the modulating frequencies are identical. 


The guiding parameter is the total excursion in phase of the 
Alfven wave which , in the language of modulation theory, is 
the phase deviation. As the modulation frequency is increased, 
the side hand structure will expand over the frequency domain 
hut the number of relevant side hands is uniquely determined 
hy the phase deviation. It is not expected that for Alfven 
waves that this will ever exceed three side hands. 

The situation is more complex when tubes of force which 
have skewed or opposed polarities sweep across the magnetom- 
eter. For, in this instance, the instrument sees a time rate 
of change whose spectrum is to he represented hy a Fourier 
integral . 

4.2 The Circ ularly Polarized Mode 

This case has no essential differences from that discussed 
in the last section except that the time variation in | now 
becomes important enough to warrant consideration of the mixing 
of this tone with the spectrum from the vector modulation. 
Earlier the angle of the plane of polarization was assumed 
nearly time invariant which was taken to include slow rotation 
of the plane of polarization. 

For simplicity it is assumed that the spacecraft is aligned 
with the preferred direction of the field; that is , one of the 
axes, say 5 3 , corresponds to the ambient field direction. 

For the wave propagating above the gyro cutoff, i is equal to 
the gyro frequency for the ions under consideration, which we 
take here to he protons . The components of the field in the 



(X1PC2PC3) set are as in eqs . (12) with £ periodic. The 

signals from the xi and $2 sensors are in quadrature as are 

— > — > — > 

the components. Hi and H 2 . The field component along H 3 is 

fixed and represented in eqs. (13) hy collapsing all the terms 

to zero except J 0 . From eqs. (13) it is seen that terms of 

the type 


and 


cos | sin(2n + l)vt 


sin £ sin(2n + l) vt 


are generated, respectively, for the Hi and H 2 signals. 

Thus | heats with an odd harmonic spectrum, generating addi- 
tional side tones. The final side hand composition is formed 

from the expansions of eqs. (13) and are 
00 

a Hi = H c ^ J 2n+ i(5) X [sin[(2n + l)v + £]t + sin[(2 n +l)v - £]t] 

n=- 00 

(l4 a) 

and 

00 

= H q ^ J 2 n+ i(S) X {cos [(2n + l)v - £]t - cos[( 2 n + l)v +£]t] 

n=-oo 

(l4h) 

The a Hi signal is composed of line pairs above and helow £, 
spaced at intervals of ±(2n + l)v while A H 2 contains side 
hands of the same form hut with phase in quadrature with the 
a Hi terms. The meaning of these expressions is straightforward. 
The amplitude of the side hands is given hy the J 2n +i terms 
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and the side hand phases oscillate in quadrature with respect 
to the field components along each of the axes. The two sensors 



have a side hand relation displaced by jt/2 since when the 
xi signal is peak; its rate of change is zero while that for 
the x 2 sensor is a maximum • The vector sum of each of the 
side band contributions from these two sensors has constant 
amplitude and contains only the angular frequency (2n + l)v. 
(See fig. 9.) 

4.3 Modulation of the Search Coil 

The apparent simplicity of the search coil as a field 
detector is subject to reconsideration in the presence of a 
hydromagnet ic radiation continuum. This instrument as orig- 
inally used measured a projected value of the ambient field 
together with time variations. The ambient field could be 
detected as a carrier tone since the coil was used with a 
rotating spacecraft. Later usage has employed attitude sta- 
bilized craft; resulting in a simpler measurement problem. 

For the spinning spacecraft; the problem of measurement 
of field variations is compounded by the combination of 
differentiation inherent in the detector together with the 
rotation of the spacecraft . Since this subject has been 
explored earlier; only the principal results will be dis- 
cussed here [ 7 ] • They have special interest since the dif- 
ferentiation process introduces an apparent phase modulation 
in the presence of a strictly scalar radiation field. The 
inverse process also takes place when the field contains only 


an ensemble of Alfven waves • 



I 
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Utilizing the coordinate systems specified earlier, the 
voltage output from a search coil spinning in a field con- 
taining time dependent terms along all resolved components 
is given by* 

e = _[H* cos cp +H2 sin 9 + w(-Hi cos 9 sin 9 + H 2 cos ©cos cp 
+ H3 sin 9 ) ]cos cot 

-[-H* cos 9 sin cp + H* cos 9 cos 9 + H3 sin 9 
- co(Hi cos 9 + H* sin 9)]sin cot ( 15 ) 


where the angles, 9 and 0 , are those defined before and H*, 
Ha, and H3 are the components of the disturbance field which 
together with H 0 form the total field 

H = (Hq + H*)xi + H2X2 + H3X3 (l6) 


in the magnetic set ■where the ambient field lies along Xi« 
The mixing of modulation is illustrated most simply if we 
assume that 9 = 9=0. Then the signal amplitude is given 

t>y 


|e(t)| 2 = (H*) 2 + (Ha) 2 + w 2 [H 2 + (H2) 2 ] - 2WH1H2 + 2 wHiH 2 

( 17 ) 


and the phase by 


i ’ = tan 1 


H 2 - wHi 
Hf + coH| 


(18) 


It is to be noted that mixed modulation appears in these 
formulas. For example, amplitude modulation is present when 

*It is assumed in the discussion of this section that the 
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turn-area product of the coil is unity. 



• _w • 

= 0 and phase modulation when H 2 = 0, the former 
corresponding to the case where amplitude variations in the 
field are absent and the latter where phase variations asso- 
ciated with Alfven waves are absent. 

The phase modulation of a pure compressional wave is 
given , from eq. (l8) , by 

* = tan-^HSMH 0 + H *) ] (19) 

A specific example of probable phase modulation due to an 
amplitude field taken from an early space probe flight [8] 
is shown in fig. 10. It is possible to assign the originating 
signal to a scalar time dependent field in this case since 
eq. ( 19 ) shows a dependence in the denominator upon the ambient 
field; t hus, for fields where the ambient value ^ R n , is suffi- 
ciently large, as in the case of the data of fig. 10, it becomes 
possible to assert that the inverse (i.e., an initial phase 
modulation) could not have been responsible for the joint modu- 
lation of the signal . 

For search coils it is possible to remove the ambiguities 
associated with differentiation by driving two coils in quadra- 
ture. This is shown as follows. From eq. ( 15 ) the two quadra- 
ture signals are given by 

e a = Hi + wh| (20a) 

and 

€ p = H* - ojHi (20b) 
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These egressions being linear differentials can be solved 
simultaneously, yielding 

Hi = (e a - W6 p )/(s 2 + ua) (21a) 

and 

Ha = (e p + we a )/(s 2 + w 2 ) (21b) 

where 

s 2 = d a /dt 2 

5*0 Multirate Switch ing and the Generalized Problem of Alias 

From the cases of modulation discussed earlier, especially 
vector modulation, it is clear that for accurate and undistorted 
reproduction of large amplitude Alfven waves the bandwidth 
extension becomes far greater than for simple scalar modulation 
of the magnetic field. Truncation of the side band structure 
in order to satisfy the requirements of the sampling theorem 
can lead to distortion from the elimination of significant side 
band energy. 

Under the simple conditions imposed by periodic sampling 
the application of the sampling theorem to band limited data is 
straightforward. When multirate switching is employed, the 
sampling problem is greatly increased in complexity. Such 
sampling is commonly found where data channels are time shared 
between experiments. That is, the signal channel is alternated 
between several data sources, with the result that each may be 
sampled at a high rate for a short time followed by a period of 
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quiet while the channel is turned to another source . The 
fundamental result of multirate sampling is to complicate the 
foldover problem greatly as will be brought out later. 

5 *1 Single Rate Sampling 

As a prelude to the mult irate case ; we now review the 
simple periodic switch. The switch function is represented 
by the sequence of deltas 

oo 

g(t) = ^ 5(t - nTi) (22) 

n=-oo 

repeated at the basic time, Ti. The Fourier transform is the 

analog of the time function in the spectral domain 

00 

G(w) = Tl X 5(w " nWs); Ws = §1 (23) 

n=-oo 

Given a sampled time series of data represented by f (t) , the 
spectrum of the sampled data is then represented by 

oo 

F( u) . £ F(io - nw g ); w s = |S (24) 

n=-oo 

which means that the original spectrum is duplicated at inter- 
vals of 2 jt/Ti and bears little resemblance to the original, 
except for the form of the replicated spectra. The function 
F(oo) is shown diagrammatically in fig. 11. Aliasing arises when 
the data spectrum is extended so that overlap takes place. In 
particular for real systems, low pass filtering is employed so 
that the overlap between the lowest two spectra is not significant. 
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although in principle "bandpass filtering about higher frequency 
replicas of the original spectrum would yield the same limita- 
tions upon folding. It is clear that the faithful reproduction 
of hydromagnetic waves must take into account the folding proc- 
ess^ and that simple band limiting to guard against overlap may 
result in side band truncation, especially when a complex PM 
spectrum is present. If the system is designed for scalar 
modulation cutoff at some frequency, v, aliasing will be elim- 
inated at the potential risk of vector modulation side band 
elimination and distortion of the final signal for Alfven waves 
of frequency ~v, unless allowed for in the design. For the more 
complex case, for example, where spin modulation is superimposed 
upon a complex PM spectrum, spectral folding can cause distor- 
tion as spectral components from n > 0 may fold into n = 0 
even though band limiting is carefully considered. 

Since the sequence, g(t), represents an idealized case 
where all the Fourier coefficients are of equal amplitude, the 
spectrum extends to infinity. In the real instance, finite 
width pulses are employed (fig. 12). Then the coefficients are 
altered according to 


C = T sin flkT/Ti 
n T l nk-r/Ti 


( 25 ) 


where t is the pulse width and Ti the repetition time, and 
the spectrum of the sampler g(t) becomes 


oo 

- 7 f B(u - 

l— j jikr/Ti 

k=-oo 


(26) 
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where oj s = 2rt/Ti, and we call the finite width sampler, g*(t) . 
There is no special significance to this result since interest 
centers on the lowest frequencies and generally the formula 
with constant coefficient suffices, 

5-2 Multirate Switching 

This problem is dealt with by one of two methods. Either 
the second switching spectrum,* which is now introduced, is 
taken as a modulation of the original singly periodic function, 
g(t) , and the spectral forms are combined, or the function 
g (t) is decomposed into separate parts whose phase relation 
is considered. In either case the essential question to which 
the discussion is addressed becomes the foldover modification 
to the low frequency domain of the spectrum. 

It is desired to find the Fourier transform of the doubly 
sampled function represented by 

tf(t) • g(t)] g *(t) = Jit) • g *(t) (27) 

where g(t) and g*(t) are defined in the previous section. 

(The switch with finite pulse width, the burst function, is 
taken for the function which interrupts the sampling given by 
g(t).) The functions f(t), g(t), and g*(t) are delineated 


In the following discussion, g(t) will occasionally be 
referred to as the sampling function, and g*(t) as the burst 


function . 


for the time domain in fig. 13 . According to the convolution 


theorem the transform of two functions equals the convolution 
of the transforms. Applying the theorem to the present case 
where the periodically sampled function, f(t) • g(t), is modu- 
lated hy the hurst function, g*(t), one has 

[f(t) * g*(t)] = F[f (t) ]*G[g*(t) ] = F(w)*G*(u)) 


pOC 

= F(u) - A)G*(a)GA (28) 

Substituting from eqs . (24) and (26) we have* 


n 00 

r— 0° — 

F(to)*G*(w) = / 

O_oo 

rjT^ ) F(w - nuig - A) 

- n=-oo — 


x 



sin jikr/T 2 

jtkr/Ts 



(29) 


With the assumption that proper convergence is observed for 
the sequences, the summation may he removed from the integral 
and the integral reduces to 



- A)s( A - 



dA = F(w 


nu) s 


kw r ) 


(30) 


where co s = 2jt/Ti, oo T = 2tc/T 2 j and Ti, T 2 , t refer, respectively, 
to the sample period for f(t), g*(t), and the pulse width for 
the switch, g*(t) . Then 

*In the multirate notation two periods, Ti and T 2 , appear 
for g(t) and g*(t), respectively; therefore, in writing eq. (27) 
using eq. (26), we must consider the time interval Ti given in 
that equation as our new T 2 . 



X f(w - nw s - kw T ) 

( 31 ) 


F(w)*G*(u>) 


oo 00 ^ ^ 

V V 1 f ' r ) Sln ^ kT / T 2 

l—i l—i Way jtler/T 2 

n =_oo j c= _ 00 v ^ 


The terms ahead of F represent a modulation of the familiar 
form given by the tabulated function sin x/x. This is now a 
doubly modulated function for which each value of the variable 
n summation takes place over all values of k. Conversely, 
one can take the opposite point of view and sum over n for 
each k. Which choice is made is a matter of convenience 
arising from the relative periods, Ti and T 2 . For the case 
at hand, T 2 » Ti and one finds that the original replication 
of f(t) at intervals of Ti now occurs for all values of k. 

It is clear that what has taken place is to replace the 
replicated spectrum shown in fig. 11b by a new infinity of 
replications. For each n, k ranges from -oo < k < oo. How- 
ever, the modulation from sin x/x becomes important and 
suppresses the lines corresponding to the higher values of k. 
For simplicity fig. l4 shows only the spectral lines corre- 
sponding to the multirate switch g(t) • g*(t) . This displays 
the multiplicity of lines which are generated. The case 
taken is for t/T 2 = 0.25^ which is for a burst function, 
g* having a duty cycle of l/4. 

Next we turn to the application of the above switching 
spectrum to a data function, f(t), which is arbitrarily band 
limited according to the singly periodic Nyquist limit given 
by the heavy line in fig. 15- Here a simple monotonic function 



was assumed for the original data spectrum. The duty cycle 
for g*(t) is 0.25 as before so that every fourth line is 
suppressed. Further , the ratio of repetition frequencies 
T 3 /T 2 =8 so that the second order (k = 5 , 6 , 7 ) for n = 0 
lies in the domain covered in the lower figure. These spec- 
tral components are not shown because they would further 
complicate the diagram. We have also not shown the inversion 
where negative frequencies are folded over the positive domain. 
Clearly, for the spectra not centered on zero, further compli- 
cation would ensue. Therefore, the figure is to be taken as 
a lower bound on the folding problem. The region of interest 
is picked by the original constraint that the original periodic 
(Ti) folding limit be observed. Attention is directed to the 
two basic types of alias which take place. The zeroth order 
folding about n = 0 is shown in the upper figure and results 
from the shifting of secondary spectra generated by k = 1,2,3. 
Since all the energy comes from multiplets of n = 0, it is 
apparent that this becomes a basic property of mult irate 
switching. Actually the same argument applies to the first 
order shifting for the multiplets of n = 1. In this case 
much of the energy resides above the intrinsic singly periodic 
Nyquist limit and does not enter into the pass band. However, 
some energy is folded down and must appear in the data band. 
Little can be done to eliminate folding from the multirate 
case unless the multiplet structure is collapsed. This can 
only take place when the repetition frequency is decreased; 



in the limit this frequency approaches the original singly 
periodic rate. The actual error introduced into the shape of 
the spectrum depends in an important manner upon the true 
spectral function of the data. Thus it becomes impossible 
in a self-consistent way to determine the spectrum from the 
aliased case. 

Some general rules can be found from the above discussion. 
Decreasing the duty cycle may be helpful since this factor 
enters in the amplitude of the multiplet structure. However , 
higher order multiplets appear within the domain of n = 0 
and the first zero of the subsidiary structure is not attained 
quickly. Thus, although the amplitudes of the multiplet 
spectra are reduced, the terms damp more slowly. The other 
variable at the disposal of the designer is the sampling ratio 
for the functions g(t) and g*(t) . As indicated above, this 
implies turning on the basic periodic switch, g(t) , less often; 
this is at the expense of information rate. In this instance 
the multiplet structure shrinks in the frequency domain and 
higher order terms of both positive and negative sign, accord- 
ing to the value of sin x/x, appear within the data band. 

Some cancellation of terms will take place, and since the 
higher order spectra are down in amplitude , they will not 
make important contributions to the final spectrum. 

The situation found for the double rate switch is mathe- 
matically similar to Fraunhofer diffraction and appeal may 
be made to physical optics for an analogy. In the earlier 



discussion of spectral folding, the higher order contributions 
to sin x/x were ignored so that the pictorial representation 
would not become forbiddingly confusing. Actually, the inten- 
sities of these maxima are not especially important to the 
folding problem as they are distributed in the following way: 
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0 


1 


12 

2 


(2/3 *) 2 


5tc 

2 


(2/5*) 2 


Jh 

2 


(2/7*) 2 


5*3 The Switching Fu nction in Vector Space 

Before turning to the simulation of the concepts discussed 
earlier, an alternative representation of the switching function 
is given here. This will provide a direct means of measuring 
the amplitude of the various terms in the switch spectrum. It 
also facilitates the view of the relative importance of the 
different contributions toward a particular spectral component. 

In the framework to be discussed in this section, the 
basic argument employed is the shifting property of the Fourier 
transform. In the language of communication practice we con- 
sider an interrupted set of sampling pulses as given in fig. 1 6 
which shows j equi-interval pulses followed by a blank period, 
corresponding to the earlier case of double sampling where the 
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periodic switch g(t) was modulated by the burst switch, 
g*(t) . The j samples are contained in the time JTi 
and the basic repetition rate is given by the period, T 2 . 

The sampling pulses are taken in the narrow limit so that 
the sequence is represented by eq. (22) . 

Now consider the case where g(t) is replaced by a 
sequence of functions each of the form g(t) but time shifted. 
Thus j = 1 represents the first term of a sequence of deltas 
repeated at intervals of T 2 . Each of the remaining pulses 
in the series j is also the first member of a new sequence. 
What we have done is merely to sum the sequence, g(t), differ- 
ently as follows: 

CO 00 

•I 

n=i n=i 



5(t - A - nT 2 ) + . . . + ^ S(t - jA - nT 2 ) 


where the repetition interval Ti is replaced by T 2 , qA is 
the time lag for the qth subsequence where 1 < q < j and 
A < T 2 ; therefore 

j 00 

&(t - qA - nT 2 ) (32) 

q=i n=i 



From the shifting theorem the companion transforms are given by 

2Jti JAjO 


g(w) = A 
v t 2 


, v 2JtiAw 

6(w - nio s )e + 


♦ . . + ^ 6(w - 


nu) s )e 


n=-oo 
J 00 


n=-oo 


“ill s( “ - ““s ’ 6 


2 itiqAw 


q=i n=-oo 


( 33 ) 



where q and A are as defined in the case of g(t) . 

Therefore each spectral line from each of the sub-sequences 
has a characteristic phase. It is the vector addition of 
the contributions from the different sub-sequences to each 
of the spectral components which determines the relative 
strength of that line . Thus for the cancellation of certain 
spectral components from the switching spectrum it becomes 
clear that the procedure to follow is to make all the contri- 
butions from the different sub-sequences cancel. 

It is possible to represent these contributions by a 
vector plot which will be unique for each of the spectral 
components. An example is given in fig. 17 which shows the 
distribution of components from the different sub-sequences. 
The problem is similar in certain respects to that encountered 
in the simple theory of diffraction. It is seen immediately 
that for cancellation of a certain spectral line, the vector 
distribution over the circle of radius l/Ti must be uni- 
formly dense. Further, the largest contribution to a spectral 
line takes place when the vectors are distributed over one 
quadrant or less . Distribution over more than one quadrant 
tends to damp the net strength of the line . Each of the har- 
monic components in the total spectrum is determined by its 
vector diagram in like fashion. Since the contributions 
become periodic in times greater than T 2 , it is necessary 
only to consider the contributions generated from the first 



members of the train of sub-sequences* In the synthesis the 
reordering of g(t) into sub-sequences means that for each 
sub-sequence the characteristic time becomes T 2 rather than 

Ti- 

In the representation given in this section, the funda- 
mental view is that one begins with a sequence having a char- 
acteristic time, Ti and removes first one pulse at a repetition 
period of T 2 . This immediately introduces spectral components 
at the frequency characteristic of T 2 > namely 2jt/T 2 . These 
components shrink or grow according to the arguments given here . 
The relative magnitude of the various contributions to the total 
spectrum can be obtained by a simulation as is shown in fig. l8 
which gives the events when a periodic sequence is interrupted 
by removal of members of the set whose basic frequency is 
5 in arbitrary units. A new periodicity equal to 1/5 is 
introduced . 

6 .0 Sequential Vector Sampling 

So far all the arguments developed have tacitly bypassed 
the fact that measurement of magnetic fields is customarily 
carried out by sampling the three components separately. Any 
proper measurement will, of course, require that both the field 
magnitude and direction be known so that three quantities are 
required in the general case, whereas the earlier discussion 
assumed that a scalar time series was involved in the data 


transformation process. 


It is assumed first that we are dealing with a time 
stationary spectrum and that aperiodic phenomena such as 
sudden changes in the field are absent. Later the transient 
case will be considered. For the purpose of this discussion 
the magnetic field is divided into three orthogonal components 
whose values individually are given as three time series linked 
by the field equations. Sampling of each of the component time 
series is taken to occur as before, using g^(t) where 
i = 1, 2, 3* The switching functions g 2 and g 3 are shifted 
in time, respectively, from gi by Ti and 2Ti« 

In the limit of periodic sampling, measures are taken at 
the rate 2 jt/Ti; however, for any one component, the frequency 
interval is 2ot/3Ti. In the other extreme, an additional 
modulation or blanking function g*(t) is added to the ensemble 
where the on-time for g* is arbitrarily taken to be 3 Ti* 

This corresponds to sampling the total field at intervals T 2 , 
each sample comprising the three components separated by inter- 
vals of Ti* (The addition of more samples during the on-time 
of g* would lend little to the discussion and is therefore 
not considered.) 

For either of the two limit cases proposed, the sampling 
is divided as in the last section into sub-sequences, each 
sub-sequence representing measurements from only one sensor. 

It is clear then that the basic sampling rate for each sub- 
sequence is given by 2at/T 2 independently of Ti and that 



the latter quantity has no bearing on the proper value for 
Nyquist band limiting. Further, if this alias-free limit is 
to be observed, then, for either extreme given above, the 
Nyquist limit categorically removes the possibility of sample 
error due to nonsimultaneity of sampling. This would clearly 
not be the case if the constraints stated before were violated. 

The correlation functions or convolution integrals are 
related to the transforms in the following manner 

r»QO 

H ± (t)Hj(T - t)at = F[H ± (t)]F[H (t)] (3^) 

where H^, Hj, and are the field components and F repre- 
sents the Fourier transform. It is not possible for the band- 
limited correlations to be affected by the switching sequence 
for their associated Fourier transforms are band limited by 
the alias-free requirement. Thus sequential sampling is 
allowed without data distortion under the guidelines laid 
down above. This in no way absolves the earlier limitations 
regarding multirate aliasing. The allowance of sequential 
sampling in the presence of multirate alias and band shifting 
requires examination but is outside the scope of this paper. 

We now turn to the case of transients in the field. It 
is simplest to take the case of the step function. We must 
define the step according to whether it is a change in the 
magnitude or direction. The step will be taken to be a contact 
surface such that the magnetic field can excite a jump in 
magnitude for the simple case where no change in direction 



(skewing) takes place. The component measures of the field 
all vary together and are linked simply hy the resolution of 
the field. Thus the response of the three detectors will lag 
in time according to eq.( 32 ) and their transforms hy the phase 
as expressed in eq. (33) • The important thing to observe is 
that the time domain is apparently reversed because the last 
detector to see the discontinuity becomes that reading the 
signal nearest the origin. Because of the phase shifts in 
the spectra, the contributions from the three sensors must 
be combined by observing the complex coefficients of the 
Fourier integral of the discontinuity . Thus, for example, 
the net amplitude spectrum is computed using the complex 
arguments. However, again as before, if the Nyquist cri- 
terion is observed, important contributions to the spectrum 
cannot be generated in the domain where phase shifting is 
important. For the low frequencies below the Nyquist cutoff, 
the phase shifts become small and the contributions from the 
three components of the field tend toward zero phase differ- 
ence . The error will diminish as the frequency cutoff is 
lowered according to eq. (33) • The allowable error in the 
spectra is then the determining factor whether sequential 
sampling is permissible. 

7.0 Synchronous Demodulatio n From Spinning Spacecraft 

In order to remove one source of modulation from data 
obtained from spinning spacecraft, it is possible to demodulate 
using an inertial or nearly inertial reference such as the sun. 



Such a procedure requires that the generation process contain 
information which is in time quadrature. A particular example 
applies to magnetometers wherein the sensors are so arranged 
that they generate signals which are phased by 90 ° , as when 
orthogonally positioned on the spacecraft. When the two sig- 
nals are each multiplied by sine and cosine replicas of the 
spin signal, the rotational modulation can be shown to cancel 
completely with a reduction of the signal spectrum to that 
obtained with inertially stable spacecraft. 

It can be shown easily that an arrangement as given in 
fig. 19 , where both the sine and cosine outputs from magnetom- 
eter sensors are each multiplied by sine and cosine terms at 
the spin frequency, will suffice to remove the spin modulation 
from the signal spectrum. Simulation can be made 
using an analog computer where the machine is fed a set of 
quadrature signals which functionally represent that obtained 
from a set of orthogonal magnetometers rotating in a static 
field. To the input are then added either or both amplitude 
and phase variations which represent the physical situation 
where compressional and/or Alfven waves are present. Such a 
simulation demonstrates that any frequencies associated with 
field fluctuations are recovered even though they may extend 
through and beyond the spin frequency of the spacecraft . The 
primary restriction becomes that of the internal bandwidth of 
the demodulation system. This must be able to transfer fre- 
quencies well beyond those representing the variations in the 
field because of the generation of a rich side band structure. 


In figs . 20 through 22 are shorn the sensor output from 
one axis, the demodulated sine and cosine outputs, and the 
associated spectra as determined from a wave analyzer. The 
case shown is scaled up in frequency for convenience . It 
represents an Alfven wave having phase excursions of ±jt/2 
radians at a frequency of 5 cps superimposed upon a spinning 
system with a rate of 1 6 cps. The appearance of these spectra 
is different than those cases hypothesized earlier where the 
negative frequency domain was employed. Thus in the present 
case the spectra are found folded according to the way the 
power is distributed. No loss of generality ensues but the 
interpretation becomes less straightforward. Figure 23 shows 
the Lissajou representation of the demodulated wave correspond- 
ing to .the vector or phase modulation of the original ambient 
field. 

Lastly, to demonstrate the recovery of signals under 
extreme conditions, an especially pathological time function 
was introduced into the computer. In this instance the spin 
frequency was half that of the modulation. Both amplitude 
and phase modulation were simultaneously superimposed upon the 
static field. The frequencies were identical for reasons of 
simplicity of recording. The amplitude variation varied the 
base field from zero to 2B Q causing extreme variations. The 
phase variations were through ±?t/2 radians. The impact of 
such extreme excursions upon the field is shown in fig. 2k- 


which is a Lissajou figure of the signals taken from the output 
of the quadrature sensor simulators. The undriven signal would 
have heen circular , representing the rotating field. Single 
axis waveform is shown in fig* 25 and the original input repre- 
senting the nonrotating field is shown in fig. 23 , which accurately 
represents the output of the spin demodulator system. Thus, 
complete recovery of this pathological signal at twice the spin 
frequency was attained. 

Demodulation of the type discussed here need only he carried 
out on the spacecraft when the maximum bandwidth is insufficient 
to transmit nontruncated signals to earth. From the earlier 
discussion, it is clear that spurious signals can, in principle, 
always be generated over a restricted range in frequency for 
any system. To avoid this it is essential that space demodu- 
lation be employed. 

8 .0 Conclusions 

The general problem of the correct reproduction of the time 
series created by signals from magnetometers which detect 
hydromagnetic radiation together with the correct specification 
of the spectra of such fields is found to place two general 
constraints upon the design of magnetometer sensor systems. 

These have to do first with the generation of modulation from 
the beating of the signals with the spin of the spacecraft in 
the event the latter is stabilized by rotation, and secondly 
by the process commonly employed where the telemetry system 
is shared in time multiplex fashion between several experiments. 
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It is obvious that simple design constraints can eliminate 


these problems. However , the conceptual simplicity of these 
constraints is not usually met in practice due to a variety 
of circumstances mostly having to do with weight and communi- 
cation limitations. To eliminate the problem of spin modula- 
tion of either type discussed in this paper requires merely 
that time quadrature signals be returned from the spacecraft 
and demodulated synchronously at the receiving station. The 
scheme shown for use in spacecraft is equally applicable on 
the ground. Generally such a practice cannot be carried out 
because of communication limitations; that is, the signal 
bandwidth is insufficient to transmit the side band structure. 
It Is clear that demodulation aboard the spacecraft is more 
economical of bandwidth in any instance. Further, the gen- 
eral requirement that quadrature signals be available is 
mandatory for demodulation either on the ground or in space, 
irrespective of the available bandwidth. (For search coils 
the requirement of quadrature sensor signals is present for 
any frequency as a result of the inherent differentiation.) 

For vector or Alfven wave modulation, the problem becomes 
more complex. A practical rule to observe is the employment 
of synchronous demodulation whenever a hydromagnet ic radiation 
field is expected to be encountered. 

In regard to the second basic problem posed, that of 
multirate sampling, here the general rule to observe is to 


periodize the data sampling procedure and employ buffer 
storage in the event that the telemetry system is time shared. 
The other recourses are either to adjust the switching functions 
for minimum tolerable folding, which is expensive of bandwidth, 
or to use multichannel telemetry such as the classical FM/PM 
system employed in early spacecraft. 

As noted early in this paper, questions generated by the 
existence of spin modulation apply equally to other forms of 
data collection where angular scanning is employed and physi- 
cally real conditions exist for the presence of data frequencies 
that encompass the spin rate . 

In conclusion it should be pointed out that the most gen- 
eral type of radiation field will be composed of a mixture of 
Alfven and magnetosonic waves. Under the assumption that the 
system response is linear, the side band structures from the 
different types of modulation can be linearly superimposed. 

Since the amplitude spectrum contains only two components which 
are phased with the carrier and these are of the same frequency 
as the first two side bands from the vector modulation, the 
spectral components merely add algebraically. The final modu- 
lation is made asymmetrical but is unfolded in the spin demodu- 
lation process. Therefore, so long as linearity is preserved, 
the mixing introduces only academic considerations. 


Lastly., although the problem of the spectral content of 
hydromagnet ic signals can be dealt with in a straightforward 
manner by calculation in the spectral domain, inversion back 
to the time domain in which measurements are carried out 
becomes difficult because of the loss of phase information 
when transforming. Phase is not recovered by the inverse 
transformation. Though detailed investigation of short 
swaths of multisampled data is possible and deductions made 
regarding wave properties, the problem of spin modulation is 
still present for spin stable spacecraft and cannot be dis- 
regarded even in the time domain. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., April 19 ., 1965 
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Figure 1.- Surface of constant phase for the fast and slow modes. The quantity A is the Alfven 
velocity, C is the sound velocity, and Q is defined as the smaller of the sound and Alfven 
velocity. The semiminor axis of the outer ellipse is given by the larger of the two velocities. 
The Alfven velocity surface is not shown. 





Figure 3.- Vector representation of the terms of eq. (7). The two modulation vectors represent- 
ing the_side bands rotate oppositely and with equal velocities with respect to the frame in 
which H 0 is stationary, at the rates ±w. Thus the quadrature components cancel always and 
the modulation is strictly amplitude. If w > w s the lower vector rotates in the backward 
direction with a higher rate than Hq. Nevertheless, the phase coherence between the two side 
band vectors is preserved. 
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Figure k.~ AM side bands for the three cases discussed in the text, the first being where 
w « u s , the second where w s < w < 2w s , and the third case where w > 2w s . Folding about 
zero frequency introduces asymmetrization of the side band structure. The lower side band is 
shown both in the negative frequency domain as well as folded. 
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Figure 5a.- View of the asymmetric truncation of the side hand structure for amplitude modula- 
tion when w » u) s . The case shown is for a system filter response with finite slope . 



Figure 5h.- Qualitative illustration of how phase modulation is introduced by the truncation 
process. The side band vectors are no longer of equal amplitude , with the result that H 0 
is phase modulated with the instantaneous angle e. 
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Figure 6.- Results of the as symmetrical side hand vectors of fig. 5 into in-phase and quadrature 
components involved in the frame of reference in which H 0 is constant. 
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Figure 8.- Side band spectrum for a monochromatic Alfven wave having the deviation discussed in 
the text . The upper figure represents the structure for Hi^ 2 and the lower for H3 as 
given in eqs . (13) . Components beyond J 3 are below 1 percent of the carrier in amplitude. 
Since the modulation is about zero frequency, the spectrum extends into the negative frequency 
domain, the lines fold over and account for the factor of 2 appearing in eqs. (13) . In 
ordinary phase modulation the spectra are distributed about a nonzero carrier. Thus the fact 
that the odd side bands are in time opposition about the carrier produces a condition differ- 
ent from the above case where the spectra fold about zero and add in power in spite of the 
opposed phase . 
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Figure 9*~ Upper diagram shows the splitting of the spectral components for Hi^ 2 from the 
beating of the spin frequency with the Alfven wave spectrum. Each line is half the former 
amplitude. The lower figure shows the spreading and asymmetrization of the spectrum for the 
case where the spin frequency is of the order of the modulation frequency. Mapping about 
zero frequency occurs as the lines shift position. Whereas in the upper spectrum the right- 
and left-hand sides are in time opposition, in the lower spectrum the phases are scrambled 
with respect to zero frequency. Thus the ensuing time function is greatly complicated as the 
vo power is not simply additive as in the case cited in fig. 8. 
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Figure 10.- Hydromagnetic spectrum of the earth* s field at 0 - 7*5 radians 
per second taken with a single search coil at approximately 7 earth 
radii on October 11, 1958 . Both the integral and differential spectra 
are shown. It is suspected for the reasons given in the text that the 
strong line at 5 radians per second in the phase spectrum (fig. 10b) 
is due to pseudo Alfven wave generation [8] . 
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Figure 11.- Upper figure denotes time domain for the function f(t) and the sampling function, 
g(t) . Lower figure shows the replicated spectra for the case of the idealized delta function 
switching case. The case shown displays some alias. 
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Figure 12.- Spectrum of the sequence consisting of an infinite train of finite width pulses 
(hurst function) . The amplitudes and phases are modulated by the term 
(t/Ti) [ ( sin 3 tkT/Ti)/(rtkT/Ti) ] . Convolving this function with the periodic switch, G(w), 
produces a sequence which gives the multiperiodic spectrum. Note that in this figure the 
pulse interval is Ti« When developing eqs. (28) through (31) the above function interval 
is labeled T g to distinguish it from the interval of the periodic switch, G(^) , whose 
interval is time Ti« 
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Figure l4.- 
extended 
t/t 2 ~ 0 


■ Multiplet generation showing only the first lobes of sin x/x for clarity. Actually 
values display phase reversal which affects line amplitude. The case shown is for 
25. The negative gj axis is not shown. It is the mirror of the above. 
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Figure 15*- Reconstruction for a hypothetical data spectrum of the folding configuration for a 
double periodic sampling function where the ratio of repetition frequencies is Ti/T 2 where 
Ti is the periodic sampling interval and T 2 the hurst rate. The duty cycle for the hurst 
modulation, g*(t) , is 0.25 so that every fourth line is zero. The second-order multiplets 

for n = 0 are not shown. They would lie in the lower figure and he superimposed upon the 

lines shown. The inversion of the negative frequencies about n = 0 are also not shown. 

The subscripts indicate the order of the spectrum; that is, n, k = 1, -3 is the third 

negative multiplet for n = 1. 
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Figure l6.- Format for a standard telemetry distribution as employed in spacecraft showing a 
hypothetical assignment of j samples followed by 1 - j blanks. Lower figure shows the 
time distribution of switching pulses. 
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Figure 17*- Vector diagram for the distribution on the circle of radius l/Tj. of contributions 
to the spectral line at the frequency w. A similar circle is required for each spectral com- 
ponent. As the order of the spectral line is increased, the distribution of vectors expands 
and becomes modulo 2n. 
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Figure 18.- Spectra of interrupted pulse trains. The top figure represents the case where the 
original periodicity at to = 5 has been changed "by the removal of every fifth pulse. This 
causes the introduction of the lower frequency lines. The following cases are for the removal 
of 2 and 3 pulses and the last is for the removal of 4 out of 5, leaving a basic periodicity 
l/5 of the original. 
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Figure 19 .- Simplified block diagram of a spin demodulator system for use with a triaxial flux 
gate sensor magnetometer. The reference oscillator is driven with a phase locked loop 
synchronized to the spacecraft sun pulse generator. 
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Figure 20.- Waveforms associated with the simulation of a spi nn i ng spacecraft in the presence 
of an Alfven wave with an angular excursion of ±jt/2. The upper figure shows the waveform 
generated on one sensor and the lower two waveforms are those from both the sensors mounted 
in the spin plane of the spacecraft after spin demodulation. 
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Figure 21.- The spectrum of the Alfven modulation superimposed upon the spin frequency. The 
spin frequency is 1 6 cps and the modulation is 5 cps. This actual spectrum is modified from 
those shown earlier as follows. The spin frequency is high compared to that of the modula- 
tion. Thus the spectrum shown in fig. 8 are folded. Then the effect of the spin is to split 
the lines as in fig. 9 . Both even and odd lines are found for this case due to the assumed 
orientation of the wave with respect to the sensors . 
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Figure 22.- The demodulated phase spectra associated with both spin plane sensors. The 
difference in amplitudes is due to the phasing of the plane of the wave as discussed for 
fig. 21. The amplitudes of the side bands are relative to the unmodulated carrier. Note 
the folding about zero frequency of the fourth harmonic line on the low frequency end. 



Figure 23.- x-y plot (Lissajou figure) of the ambient field with vector and scalar modulation 
as it would appear in space or as measured by a nonrotating spacecraft. This is also the 
output of the demodulator. The multivalued trace is due to small errors in the function 
generators . 
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Figure 2k . - Lissajou figure of the outputs from two orthogonal sensor simulators for the case 
of simultaneous Alfven (vector) and amplitude (scalar) modulation of the base field at twice 
the spin frequency. The signals shown display the modulation superimposed upon the spin 
frequency of the spacecraft. 
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Figure 25* - Waveform of one sensor shown to display the pathological character of the net 

function. 
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